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In the framework of Standard-like maps, we find an efficient scalar algorithm to compute finite 
time Lyapunov exponents (FTLE) and Oseledets' splitting (or covariant Lyapunov vectors, CLV); 
by connecting such procedure to invariant manifolds for generic non- fixed points (covariant curves), 
we find explicit approximations for tangent vectors (i.e. CLV) and manifold curvatures. In the 
appendix we generalize the numerical algorithm to any different iable map of the plane. 

I. INTRODUCTION 



In the huge literature concerning the study of bidimen- 
sional maps, much effort has been devoted to the study 
of Lyapunov exponents and invariant manifolds, the for- 
mer being the primary indicator for chaos, the latter be- 
ing the actual skeleton of phase space structure; negative 
and positive Lyapunov exponents compose the key in- 
gredients to observe stretch-and-fold mechanism in maps 
of the plane, quantify KS-entropies in both time direc- 
tions and give information about phase space uniformity, 
just to name a few. Together, stable and unstable in- 
variant manifolds are far more subtle objects, encoding 
the geometry of motion; indeed, they represent the key 
ingriedents for a multitude of fundamental issues such as 
invariant measures and symbolic coding of the dynamics 
(construction of generating partitions). 

In this paper we first analyze the informations we can 
get upon the skeleton of the dynamics in a generic point 
of the phase space of a Standrd-like map: the relevant 
observation is that to any point we can associate two 
covariant curves (that are not generally invariant unless 
the point is periodic), whose directions are parallel to 
the covariant Lyapunov vectors (CLV) [i, Q, that deter- 
mine the corresponding Oseledets^ splitting jsl, 1^. Then 
we consider how the splitting angle (between CLV) can 
be studied in detail for such maps: we comment on its 
relationships with large deviations for finite time Lya- 
punov exponents (FTLE) [5|49|, and elaborate upon an- 
alytic approximation schemes, which are shown to yield 
fairly accurate pictures. In the final part of the paper 
we generalize the aformentioned algorithm to calculate 
CLV and curvatures of covariant curves to any differen- 
t iable map of the plane; althought for very complicated 
maps the method can become slightly cumbersome in the 
calculation of curvatures, it remains scalar and compact. 




Figure 1. Image of 10^ points equispaced along the torus 
diagonal under 8 iterations of the Chirikov- Taylor Standard 
map, as in eq. ([T]) with f{x) — 2x -\- i^sin(x) and K — t^. 



II. STANDARD-LIKE MAPS 

A. Choice of representation 

We consider the class of symplectic standard-like maps 
defined on the torus, V x = (x, ?/) G T^ : 



f(^,^) = (/(^) -y I 



modP 



, x) 



(1) 



With f{x) e &{R) s.t. /(x + P) U^p = fix). Such 
maps are always reversible [7]; namely there exists an 
involution X (a map with X^ = id) s.t.: 



Xof 



f-^oX ; 



(2) 
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notice that here: f ^{x,y) = {y , f{y) - x l^^dp)- 
Our choice of a two coordinates representation yields a 
particularly simple form of involution: X(x,?/) = (^,x), 
linear with matrix form X = [ ? q ] • 

Reversibility by reflection implies that fixed-points al- 
ways lie along the diagonal; moreover, if x* is a pe- 
riodic point also Xx* is periodic with the same pe- 



riod. This can be seen as a choice of a two coordi- 
nates representation of the ordinary standard mapping: 
(x,p) -^ {x -\- p -\- G{x)^p + G{x) , so the function 
G{x) = f{x) — 2x can be physicahy interpreted as an 
external conservative force acting on a point mass whose 
position is bounded on a circle; fixed-points then corre- 
spond to equilibria of the system, where the force van- 
ishes: G(Xfix) = , 7/fix = ^fix- 



B. Stability 

The Jacobian matrix of map ([1]) naturally reflects the 
reversibility property: 



9x 



f'{x) -1 

1 



(3) 



where the prime symbol stands for derivative w.r.t. x, 
while: 



j-^ 



-A.J X-^ 



(4) 



We may indeed flnd a more general symmetry, concerning 
the stability along trajectory segments. First we deflne 
Fby: 



dx 



FJx) 



J f fc 



Hx) 



Jf(x)Jx 



and notice that it enjoys the cocycle property: 

F,+fe(x) = F,(f'=(x))Ffc(x) . 
Since Fo(x) = 1, we see that: 

F_,(x) = F,(f-^(x))"' = 

- *Jf-fc(x) • • •*^f-i(x) ' 

and by employing (|2]) and (|5]) we may rewrite it as: 

F-/c(x) = XJ^j.fc^^^^ . . . Jxf(Xx)X ; 



(5) 



(6) 



(7) 



(8) 



by observing that, according to our two coordinates rep- 
resentation: 



^Xf^Cx) 



J f fc 



Ux) 



(9) 



we finally get: 



Ffc(x) = XF_fc(Xx)X , 



i.e. the Jacobian matrix of any map iterate in x and its 
corresponding inverse in Xx are connected by the linear 
transformation X. 

Since map ([!]) is area preserving, det(Ja;) = 1 Vx, sta- 
bility of flxed point x* depends only on the trace of the 
Jacobian matrix, /^ = tr(Ja^*) = f {x*) and gives the 
following form for the eigenvalues x±- 



g±arccosh|/^/2| ^ |/V2| > 1 
X± = Cr* • <; ^^. arccos|/;;/2| ^ 1/^21 < 1 



(11) 



and eigenvectors w± ex [^i] , with a* = sign(/^) and 
X- — l/x+ 5 i^ th^ unstable case (|/^/2| > 1) the eigen- 
values are real and it follows from (pT]) that they cor- 
respond to the cotangents of the eigenvectors polar co- 
ordinate, w± = [sfn] (q^^) , SO that: «+ = ^ — a_ , 
the eigenvectors are symmetric w.r.t. the diagonal. The 
angle between the two vectors is called splitting angle^ 
6 = a- — a-^ , which can be written as 6> = 2(^ — a+) 
; the parabolic case \fl/2\ = 1 = x± corresponds to the 
degenerate limit: a± = j , ^ = , in which both eigen- 
vectors align to the diagonal. 



C. Covariant structures &z FTLE 

In the following a similar notation will be used for non- 
flxed points in phase space Q , to which the invariant 
manifold deflnition can be generalized, [lO|, by obtaining 
sets that in general are not invariant: 



T±(x) 



lim 



|f'=(x)-f'=(y)||=0 



but fulfill the covariance condition: 



f'=(a'±(x)) = a'±(f'=(x)) 



(12) 



(13) 



i.e. are transported along with the dynamics and, in 
the limits /c ^ ±oo, converge to the single orbit gener- 
ated by point x; under fairly general assumptions one can 
show [10| that each set IP^ can be represented by a dif- 
ferentiable curve x{s) (actually G C^(M), as the map f is) 
parametrized by its arc length s and such that x(0) = x. 
It is important to notice that property (|2]) applied to ([12]) 
allows to connect sets of the two different types: 



XJ^±(Xx) = J^^(x) 



(14) 



namely, the + set for one point corresponds to the re- 
flected — set for the reflected point. This fact suggests 
that one can infer properties of the + sets from results 
obtained upon — sets and viceversa; in the next section 
we will show a method by which this can be achieved. On 
such grounds, to keep formal expressions at a minimum, 
in the following we will consider a single set ^(x). Now we 
exploit covariance ([13]) once manifolds parametrization is 
considered. Fix a reference point xq : the family of curves 
(l^J {xn(s)}nGZ coiucidc with {J^ (f^(xo))}^^^ . Covariance 



implies that for any S[o] G M , f (xo(s[o])) lies in xi(s) , 
so there exists 5[i] G M such that xi(s[i]) = f (xo(s[o])) . 
This observation leads to deflne the sequence of mappings 
(p[n] : R ^ R such that S[n+i] = (p[n]{s[n]) and: 

Xn+l (9^[n](5[n])) = f(Xn(5[n])) • (15) 

Differentiating w.r.t. sr^i yields: 



Xn+l(5[n+l])^[n](^[n]) 



^islr.^)^n{S[n]) , (16) 



the mapping between unitary tangent vectors of the co- 
variant sets. Writing (p!6|) more compactly: x^+i(/9[^] = 
JnXn , the mapping derivative: 



^[n]{S[n]) 



dsu 



^+1] 



dsu 



(17) 



represents the local expanding factor; iterating (p!6|) k 
times from n = and defining F/e^o = F/i;(xo(so)): 



k-i 

Ffe,oXo = X/e JJ (P[q] 

g=0 



(18) 



provides a consistent way to express the FTLE via 
Birkhoff averages of the local expanding factors: 



k-l 



Xk,0 



= iin 



|F/e,oXo|| = iXl^^l^wl 



(19) 



q=0 



Generally FTLE can be calculated by using any random 
unitary vectors instead of xq: in the limit k -^ ±oo 
we then get with probability 1 the biggest/smallest Lya- 
punov exponent of the system respectively (for Standard- 
like maps they have same value and opposite signs); the 
covariant Lyapunov vectors are defined to be the only 
vectors yielding the same Lyapunov exponent under both 
temporal limits. Practically speaking, this means that 
the first and second CLV are the only possible tangent 
directions that, for unstable chaotic orbits, converge ex- 
ponentially to zero in the respective limits k -^ ^oo. 
This definition, combined with those for covariant sets 
in ([12]), prove that vectors x^ , tangent to the covariant 
sets y^, actually coincide with the Oseledets' splitting 
[sl, 3] in any point in which they can be defined. 



III. DIMENSIONAL REDUCTIONS 

A. Scalar tangent evolution 

The approach we propose here, already used in [11] 
for different purposes, exploits the isomorphism between 
special linear transformations (the 5^(2, R) group) and 
rational functions (the Mobius group, see appendix A). 
Being tangent vectors x unitary, they can be uniquely 
represented by polar angles a G [— tt, tt] ; defining the two 
coordinates covering the circle: ip = cot (a) , a = sign(a) 



COs(Qf) 

sin(Qf) 



a 



^/TT^ 



(20) 



Such coordinates can be regarded as functions of the 
arc length 5, as x is, and ultimately as phase space 
functions [12]; for this reason we will refer to function 
ip{'Kn{sn)) through the lighter notation '^n, and the 
equivalent convention will hold for similar functions. By 



inserting the Jacobian matrix (J3j) into (p!6]) . we get, after 
some algebra, the evolutions: 






(21) 



which allows us to make the local expanding factors ex- 
plicit: 



^[n] 



Vn+l 



V^n 



Vn 



Vn+l 



\^n 






(22) 



The modulus on ijjn in the last equality and the second 
evolution in (|2T]) are both due to the fact that (p can- 
not change in sign, since function ip is conjugated to f 
and thus invert ible and monotone. It is remarkable to 
notice how the one-dimensional geometric evolution, rep- 
resented by the first of (|2T]) . determines the exponential 
instability as well, as shown by ([22]) : such a relation al- 
lows us to rewrite the expression ([19]) for the FTLE: 



Xfe,o 



k—l 
q=0 



iln 



1 + ^0 



A/c,o ; (23) 



^0 



in which the approximation is justified by the fact that 
the geometric factors ?/;, althought unbounded in prin- 
ciple, cannot encode in themselves any exponential be- 
haviour in /c, so we typically expect that the neglected 
term chaotically fiuctuates to zero as ^ . Due to this fact 
the approximation becomes an equality in the infinite 
limit, in which A depends onl y o n the ergodic component 
to which the orbit belongs [10|; notice that every nu- 
merical method of calculation of Lyapunov exponents is 
biased by terms falling to zero as ^ [LYA] , so the approx- 
imation actually cleans out part of them[13|- This can 
be rephrased through an asymptotic expanding factor: 

k-l 
r/c,m = W "^m+q , \ log \Tk^ni\ = Afe,m - A , (24) 

whenever k is large enough; F fulfills the multiplicative 
cocycle property, as in ([6]), by definition: 



k-\-j,q — -'- j,q-\-k^ k,q -> 
ro,g = 1 , Ti^q = ^jjq , 



(25) 



while for A the additive version of such property holds: 

(k + j)\k+j,q = j\,q+k + kXk^q , (26) 

giving the way to evolve A/e,o through Ai^^ = log \il)q\ : 

(A: + l)A/e+i,o = Ai,fc + A:Afe,o • (27) 

We thus conclude that, to calculate one of the CLV and 
the associated FTLE, only the first of equations ([2T]) is 
needed (i.e. the a evolution can be avoided). 



B. Scalar curvature evolution 

We showed the correspondence between hnearized co- 
variant curves and the Oseledets' sphtting, it should be 
interesting to go to the second order i.e. study the curves 
curvature also; as far as we know, such problem has only 
be considered in pjj and [14], where a similar approach 
is devised for the Henon dissipative map. To compute 
curvatures it is sufficent to use the second derivative of 
the curve w.r.t. the arc length exploiting the unitarity of 
tangent vectors: 



X = QfYx 



(28) 



with g smooth and invert ible, inducing the dynamics 
along the two curves g^(^) = ^ ^^^ , meaning that: 

yt — 9^^{y) 5 i-^- fo^ points on the curves the evolution 
is scalar through g^^] this can be seen by the fact that 
the conjugation holds: 



f"°g*(y) = g°5*"(2/) 



(36) 



We think that further understanding, from functional 
theory, is needed by this aspect of standard-like maps; 
in this work we limit the choice to the paradigmatic 
Chirikov- Taylor standard map, f{x) = 2x -\- K sm{x) , 
as a numerical benchmark. 



with Y the ^ rotation matrix; inserting in the definition 
of curvature Hi gives then the relation: 



XXX 



(29) 



Now we want to provide an iterative procedure for curva- 
tures as well: to this purpose we introduce the auxiliary 
quantity r] defined as 



V 



t 

y 



= -aa{l + V^)^ 



(30) 



The last equality shows how 77 is connected to a, the 
signed curvature. If now we derive the geometric map- 
ping dH]): 

^n+l^[n] = f"{Xn)Xn + "y^n , (31) 



and then use 
obtain: 



, by dividing (J3T1) by ?/nV^n = x^, we 



^n+l = f"{Xn) + 



^l 



Vn 



(32) 



Notice that in order to recover the signed curvature a one 
would need the logical function a also (and thus imple- 
ment also the second evolution in eq. ([2T]) ): this can be 
avoided if the analysis is restricted to probe the flatness 
of covariant curves, for which logarithms of curvatures 
are better suited due to the very wide values range. 
It is interesting to point out that ip and r] can be inter- 
preted also as follows: 



, X dx 

^ = - ^ 

y CLy 



ip dip 

y ^y 



d\ 



d^ ■' (^^) 



this implies that, if a functional local representation of 
covariant curves exists in the form g{y) = {g{y),y), then: 



^ 



V 



(34) 



By restricting to invariant curves, it is possible to con- 
strain such representation to be global just for chosen 
curves (inverse problem) by the use of the invertible func- 
tional class: 



fix) = g{x) + g-\x) 



(35) 



C. Evolutions convergence 

We finally consider stability issues related to schemes 
(f2T]) and ([32]) . by perturbing ip and r] at first order and 
evaluating /', f" over a fixed reference orbit; this yields: 






-3r]nlpn^ IpJ 



Stpr> 



(37) 



triangular Jacobian matrices imply that diagonal entries 
are eigenvalues and thus diagonal entries of products of 
Jacobian matrices are products of Jacobian eigenvalues; 
this, along with (|24|) . shows that perturbations Stp and ^77 
vanish exponentially fast, with corresponding exponent 
— 2A and — 3A , and A the orbit positive exponent; after a 
number N of iterations, dependent on the magnitude of 
A, the sequences {ipn} and {rjn} can be considered as the 
actual values of ip and r] for the covariant curve on the 
semi-orbit {xn}n>Ar. We remark that all these features 
are completely equivalent to the conver gen ce properties 
of any standard matricial method, e.g. [15], [16], [l], be- 
ing solely due to the hypotesis that the reference orbit has 
distinct Lyapuov exponents. As showed in Appendix A, 
the method we propose for CLV converges exponentially 
fast, with the Lyapunov exponents difference as exponent 
(and thus it works for non- hyperbolic systems also), as 
the multidimensional method recently introduced by [l|. 



D. Analytic approximations 

Evolution laws (|2T]) , (|32]) can be also examined for for- 
mal solutions that, if not in closed form, can give at least 
some insight on connections with both FTLE and the 
function /; it is possible to write them as pseudo-affine 
transformations : 






r 
f" 



(^n) + 











(38) 



with coefficients dependence from /', f" and ip ; we de- 
fine the auxiliary matrices (see(| 



■ hk 



(-i)^rf 



i,fc 




F^ 



(39) 
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Figure 2. (Color online) Chirikov- Taylor map, K — 10; in 
panel (a), evolution under ([2T)) of the log- variance of an initial 
uniformly random ensemble V^o, over a single reference orbit 
of P = 10^ iterations, V^^ (orange) and -0" (blue). Dashed 
lines are linear fits for the variance decay exponents, both 
around 3.6, i.e. twice the FTLE value over the same orbit as 
shown in panel (b) for A^(xo) (orange/blue continuous line) 
and A_^(xp) (orange/blue dashed line). 



Inserting them in (|38]) and making use of property (|25]) we 
obtain the two solutions: 



n,0 



Z^ n-9,q 



9=1 






{y,) ; (40) 



the idea now is to move the orbit initial condition from 
xo to X-^: by (|M|) the first term on the r.h.s. of (|^Q|) 
vanishes exponentially for n ^ oo, loosing memory of the 
infinite past and thus, after a sum index change, we can 
consider the limit solutions: 



lim 



/_^ q,-q 



n,-n q^Q 






(y-,) , (41) 



to check that they fulfill the respective evolutions (|2T]) 
and (|32]) . This type of solutions is implicit since, as ap- 
parent by the definition of F, it requires the evaluation 
of function tp itself on the negative orbit of x; that's why 
we consider it only formally. Rewriting separately (HU), 
through (|24|) again, we get to a more elegant structure: 



V-lx) = ^(-r 



q -2qXq-q H 



f\y- 



q=0 



9=0 



(42) 
(43) 



If, by an opposite approach, we recursively reinsert evo- 
lution (|2T]) in itself, we can rewrite series (|^2]) : 



V'(x) = f'{y) 



1 



(44) 



f'iy-i) 



f'iy-2) 



f'iy-3) 



in the equivalent form of a continued fraction with real 
entries f and alternate signs; in this way '^(x) depends 
on the backward orbit of x and is no more in implicit 
form. Truncations of continued fraction (|^1|) corresponds 
to finite sums of (j42j); on one side, this clarifies tjj phase 
space dependence, on the other it shows its complexity: 
series (j42j) is weighted by exponentials in —2q whose ex- 
ponents are FTLE, Xq,-q^ of finite time q; this means that 
the stronger the exponential reduction is, the smaller are 
the FTLE fluctuations around A. While series truncation 
may sounds reasonable, the only source of errors comes 
directly from the FTLE deviations: if, for some very high 
(7, deviations bring the FTLE to 0, the associated f'q/fq- 
dependent term is weighted nearly as the first one. By 
the equivalence of expressions (j42j) and (|44|) , we conclude 
that the heuristic accuracy of considering truncations of 
(|42]) will depend pointwise on x via the FTLE deviations 
upon each x orbit. The lowest order approximations are: 






f\y) 



nf{y)-x) 



(45) 



^ii{x,y) = f{y) 



1 



fiV-i) 



1 



f{y-2) 



with y-i = f{y) — x and y-2 = f{y-i) — y diS can be 
obtained from the inverse map ([2|); here the |modP have 
been dropped since /' already has the correct periodicity. 
In the next sections we will show that these first three 
approximations already provide a rather detailed descrip- 
tion of subtle phase space issues: from a statistical point 
of view it will be shown that, in sufficiently unstable 
cases, ipi^ and sometimes even ipo^ is already enough to 
account for the main properties of interesting numerical 
objects as the probability densities for Ai^n = In iV^nl (i-e. 
the object over which, as we have shown, Birkhoff aver- 
ages yield FTLE) and the splitting angles On between the 
positive and negative covariant curves. 



E. Splitting angles 

As we previously mentioned. Standard-like maps writ- 
ten in coordinates ([!]) enjoy the property of reversibility 
by reflection, that in turns connects positive and negative 
covariant curves (see (p!4|) ): functionally, this implies: 

^^{x,y)= ,_} . , (46) 
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Figure 3. Chirikov- Taylor map, panel (a): K = 27r, panel (b): K = 47r; in grayscale, density plots of planes (In IV^+I, In |(i+|), 
(^,ln|d+|) and (Inl-^+l,^), obtained from a single 10^ iterations orbit. Notice that, although clear trends are apparent, 
deviations are structured and systematic as confirmed by histograms in Fig. [HandEl angle 9 is expressed in degrees. 



with ip^ corresponding to the covariant sets J^^; by using 
= a~ — a~^ and some trigonometry we have thus: 



cot{0{x,y)) 



ilj+{x,y)ilj+{y,x) - 1 



(47) 



so that it becomes apparent that ^ is a symmetric func- 
tion of phase space and can be expressed through the 
positive slope ?/^+ only. Combining the evolutions of tp^ 
in (|2T]) it is possible to obtain the relation: 



^i^n 



sin(6>n) 



sin(l9n+i) 

which, by making use of (fT9|) , leads to the following: 

sin(6>o) 



(48) 



-^^,0 + \,o 



iiog 



sin(6>/e) 







(49) 



that is, the exponents converge to opposite values as ^, 
A+ = — A~ , but, at finite times, deviations from such 
equality are induced by nonuniform hyberbolicity. Notice 
that eq. (|48|) is a natural area preserving condition. 



IV. NUMERICAL TESTS 

We numerically apply our approach to the paradigm 
of all Standard-like maps known as the Chirikov-Taylor 
map (CTM) [13, [l8|; in our notation this corresponds to 
the choice f{x) = 2x -\- i^sin(x), with K interpreted as 
the strength of the external force, playing the role of non- 
linear parameter. Our numerical analysis is organized 
as follows: we first focus on a single type of covariant 
curve, i.e. 1P+ the positive one, by examining the qualita- 
tive relationships between the effective local exponent [19j 



ln|?/;+|, the curvature a+ and the transversal angle ^, 
evaluating them in each point of unstable orbits; then we 
consider the approximations in (|45]) and compare their 
functional evaluation over a square grid with the numer- 
ical results from evolution ([2T]) , averaged over the same 
grid. With tjj evaluated from such data, we finally com- 
pare the splitting angles obtained through (|47|) . To rely 
on the numerical information extracted from evolutions 
([2T]) and (|32|) it is necessary to discard the initial con- 
vergence transient; as previously shown, this depends 
on the Lyapunov exponent of the considered orbit and 
for CTM tipically consists of a few hundreds iterations. 
Since splitting angles are calculated from both the stable 
{^p~) and unstable ('0^) CLV and to calculate tlj~ one 
has to run backward in time eq. (|2T]) , for both initial 
and final segments of the orbit must be dropped (see Fig. 
[2|). Quantities like ?/^, r] and are given in terms of com- 
plicated expressions: to get an intuitive grasp at their 
mutual correlations we plot 2D histograms by cosider- 
ing the three possible planes in Fig.s O by inspection of 
such figures it is possible to notice that the curvature 
is strongly bounded by K and by the function K/tJj^] 
this can be understood by formula (|^ and the fact that 
max|/''| = K. Moreover there are clear statistical trends 
connecting the three involved quantities: to higher effec- 
tive local exponent In I'^+l there corresponds a very low 
curvature (in Fig. [3l curvature is in logarithmic scale); on 
the other hand, very small splitting angles correspond to 
higher curvature values. These observations lead to the 
euristic conclusion that, on average, the more one covari- 
ant curve is flat, the more local instability is strong and 
splitting angle is away from zero; this is rather intuitive, 
since in regions with flat covariant curves the system lo- 
cally resembles a uniformly hyperbolic system, in which 




(see panel (b) from Fig. |4]). The statistics of local ef- 
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Figure 4. (Color online) Chirikov- Taylor map, K — 27r: panel 
(a): PDF of log curvatures In |d^| along with the same prob- 
abilities conditioned by expansive (ln|'0^| > 0, orange) and 
contracting (In 1-0^ | < 0, blue) local effective exponent; panel 
(b): the same for the splitting angle 0^ expressed in rad. Both 
quantities are calculated upon the same 10^ iterations orbit. 



splitting angles are strictly bounded away from zero. We 
remark that these considerations are purely qualitative 
and have a definitive meaning only in the statistical sense. 
This can be seen also in both panels (a) of Fig.s [5] and 
m where histograms of In I'^+l and In \cx^\ confirm that 
deviations are present and decay algebraically for both 
quantities; in particular, curvature can take values very 
much larger than K and, by computing its conditioned 
probabilities w.r.t. the sign of Inl'^+I, it is apparent 
that small curvatures correspond almost completely to 
positive local effective exponents (see panel (a) of Fig. 
[4]). Applying the same conditioned analysis to the split- 
ting angle probabilities brings to similar conclusions, the 
fraction of negative local effective exponents corresponds 
only to angles around zero, again, in a statistical sense 



Figure 5. (Color online) Chirikov- Taylor map, X = 27r, panel 
(a) : PDF of the local effective exponent ln|'0^| calculated 
from a single 10^ iterations orbit; panel (6) : blowup of the 
same graph near its maxima. Vertical lines correspond to the 

values In 1 2 lb iT I (blue dot-dashed line) and In 2 zb iT — 2±' k 
(orange continuous line), i.e. eq.s (|45|) evaluated in points 
(xo,^o) s.t. cos(^o) = ±1 and cos(?/-i) = ±^1. Correspon- 
dence between such points and some of the maxima can be 
understood from the fact that, by the hypothesis of uniformly 
randomly distributed (xo,^o), the cosine values are denser 
near ±1. Magenta dashed line corresponds to ln|i^|, the 
roughest approximation of the PDF average, i.e. the FTLE. 

fective exponents can be characterized directly by the 
approximations from (j45j): in Fig. [5], panel (b), is shown 
a blowup of the PDF from panel (a) in which we com- 
pare the positions of maxima of the distribution with the 
V^ values corresponding to the first two approximations 
in (j45j) evaluated in points such that cos(?/o) = ±1 and 
cos(?/_i) = ±''1; the coincidence between such values and 
the maxima of the graph can be explained through the 
accumulation of values of the cosine function evaluated 
over uniformly randomly distributed points (x,?/), due 
to the invariance of the Lebesgue measure w.r.t. CTM. 
Finally we check the effectiveness of the approximations 
obtained in eq. (|^ . To this end we consider a square 




Figure 6. (Color online) Chirikov- Taylor map, K — 27r: PDF 
comparison for splitting angles calculated by the four meth- 
ods as in Fig. [71 Ogrid (bold black), On (bold red), Oi (thin 
orange), Oq (thin blue). Noise depends on grid size. 



grid partition of phase space: while the numerical quan- 
tities are averaged inside each cell that contains the orbit 
point to which they correspond, the approximating func- 
tions are evaluated in the centers of each of those cells. 
In the presence of resonant islands (unaccessible by the 
unstable orbits we consider), the cells covered by them 
are considered as empty and plotted in white: From Figs. 
Hit can be seen that, while t/^o can account only for the 
gross features of phase space structures, the agreement 
between the numerical results, ijji and ?/^// is quite good 
ior K = TT and extremely good for K = 2it. This can 
be appreciated in Fig. [6] also, where a comparison of the 
splitting angle statistics is shown for the four distinct cal- 
culations of panel (6) in Fig. [71 : it is apparent that even 
the lowest order approximation, ^jjq , is enough to figure 
out the histogram structure. This result comes with no 
surprise since, as we already mentioned, the goodness of 
the approximations is better for highly unstable systems 
(i.e. when the Lyapunov exponent is greater, and for this 
system it is known to scale , on average, as In |i^|). 



V. CONCLUSIONS 

By exploiting the definition of covariant curves ^] and 
their connection to the Oseledets' splitting, or CLV, we 
devise an optimized scalar algorithm to calculate them, 
along with curvatures of such curves, for any & map 
of the plane; the structure of Standard-like maps makes 
the framework particularly simple and it allows to use 
the proposed algorithm as a functional tool to construct 
CLV in approximate form expression and to show the 
connections between geometric and metric properties of 
such systems. To give a glimpse of what can be ana- 
lyzed in such a framework, we apply it numerically to 
the Chirikov- Taylor map, showing both the possible in- 



formation coded by covariant curves(e.g. the connections 
between curvatures, local instability and transversal an- 
gles) and how the unstable cases are carefully described 
by the approximations we introduced. 
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Appendix: Mobius evolution 

To extend the evolutions (|2T]) and (|32]) to any & map 
of the plane we make use of the homomorphism between 
the 5'L(2, R) group (in which any invertible Jacobian ma- 
trix is included by normalization) and the Mobius group 
of linear fractional transformation; rewriting ([16]): 



^n+l'^n — *^ri^n -> *^n 



A B 
c D 



, X72 CX 



^rz 



(A.l) 



we define 7^ = Cn'^n ^ D^; in this way we get the fol- 
lowing transformations: 



V^. 



^jA^n 



Ani^n + Bn 



(Jn+l = Sign(7n)crn , (A.2) 

SO that the local expanding factor and FTLE become: 

fc-i 



^n = 7n- 



Vn 



Vn+l 



AM = i^ln|7,| , (A.3) 



By deriving w.r.t. the arclength we obtain the extension 
of (|32]) for ?7 = ^ to the general case: 

r]n+l = — {det{J n)r]n + {an + Ki^n + Cn^l^D) 

a„ = v„ • (DVB - BVD)„ , 

6„ = v„ • {DVA - AVD + CVB - BVC)n , 

c„ = v„ • {CVA - AVC)n , 

v„ = ['^i"], (A.4) 

with V the gradient operator and " •" the scalar product. 
The L yap unov exponents for evolutions (|A.2p , ()A.4p are 
found [20| by perturbing (?/;, rj) while keeping /', f fixed; 
by combining ()A.3p with: 

lim In I det(F/e,o)| = (A+ + A") , (A.5) 

/e^-±oo 
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Figure 7. (Color online) Chirikov- Taylor map, K = 27r; in color scale, the splitting angle (in degrees). Panel (a), (b), (c): 
evaluated using -00, V^i and V^/j as in (|45)) respectively, over a square grid of 10^ x 10^ cells; panel (d): numerically calculated 
from a single 10^ iterations orbit and averaged over each cell of the same grid as before. 



we get the exponents for (?/;, r]) as functions of A+ > A 



A^ = lim - In 

k^±oo k 



A^ = hm - In 

k^±oo k 



S^ipk 



Sr]k 



Sr]o 



-(A+-A-) , 
-(A+-A-)tA± , (A.6) 



from which we deduce that the convergence of t/j only 
requires the that orbit exponents are non-degenerate, 
A~ < A+ while for r] must also hold that: A~ < A+/2 , 
a sort of hyperbolic bound in case of purely expansive or 
contractive systems. The Standard-like case corresponds 
to: an = bn = 0, c{x,y) = f"{x) while 7n = V^n and 
A^ = -2A+ , A^ = -3A+ . 
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